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ABSTRACT 

Acoustic  oscillations  were  induced  in  a  sub-scale  liquid  rocket  engine  that  burned  liquid 
oxygen  and  gaseous  hydrogen  as  propellants.  The  oscillations  in  the  chamber  were  forced  by  a 
rotating  gear  just  downstream  of  the  nozzle  throat.  High  frequency  data  was  acquired  for  pressure 
and  velocity  via  a  pressure  transducer  and  a  magnetic  flowmeter.  The  magnetic  flowmeter  obtains 
the  acoustic  gas  velocity  by  measuring  the  voltage  induced  by  the  ionized  combustion  products 
moving  through  an  externally  imposed  magnetic  field.  A  cross  correlation  was  performed  on  the 
velocity  and  pressure  signals  to  determine  the  amplitude  and  phase  difference  of  the  two  signals.  A 
linearized  one-dimensional  acoustic  model  was  developed  to  simulate  the  mean  and  unsteady  flow 
within  the  chamber  with  mass  and  energy  addition.  The  phase  difference  between  unsteady 
pressure  and  velocity  was  determined  from  the  model  and  fit  to  match  the  phase  difference 
measured  by  the  experiments.  The  points  where  the  modeled  and  experimental  phase  differences 
agreed  determined  the  real  part  of  the  propellant  evaporation  and  combustion  pressure-  and 
velocity-coupled  response  functions. 

NOMRNCIATURE 

a  speed  of  sound  (m/s) 

B  magnetic  flux  density  (T) 

E  electric  field  strength  (V/m) 

A  frequency  (Hz) 

/-  distance  between  electrodes  in  the  magnetic  flowmeter  (m) 

mi)  mass  flux  due  to  propellant  mass  and  energy  addition  (kg/m^s) 
tip  number  of  data  points  in  average 

p  pressure  (Pa) 

pressure-coupled  response  function 
^  velocity-coupled  response  function 

(/)  power  spectra 
T  temperature  (K) 

Tf  flame  temperature  (K) 

^  voltage  induced  across  electrodes  in  the  magnetic  flowmeter  (V) 

U,u  velocity  of  combustion  gases  (m/s) 

X  distance  along  length  of  chamber  (m) 


Greek  Symbols 

a  proportionality  constant  for  the  magnetic  flowmeter 

A  r  non-i  sentropic  temperature  change  (K) 

A  ^  phase  error  (°) 


V  r‘ohf»rpnf'p 


OTIC  QUALE?  m 


1 9971 1 1 8  043 
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Superscripts 

mean  property 

unsteady  property  as  a  function  of  x 

Subscripts 
p  pressure 

pu  pressure  and  velocity  correlation 
u  velocity 


INTRODUCTION 


Prior  research  has  proven  the  validity  of  several  techniques  for  analyzing  combustion 
instabilities  in  liquid  propellant  rocket  engines  including  the  time  lag  concept  and  models  of  the 
actual  combustion  mechanisms*.  Yet  another  technique  defines  response  functions  to  determine 
combustion  stability  within  the  rocket  chamber.  When  analyzing  the  stability  of  rocket  engines, 
researchers  determine  the  pressure-  and  velocity-coupled  response  functions.  Each  response 
function  is  a  complex  number  which  determines  the  response  of  the  various  processes  within  the 
rocket  engine  with  respect  to  the  unsteady  pressure  or  velocity. 

One  method  of  determining  the  response  functions  is  to  measure  the  unsteady  properties  of 
the  gases  inside  the  rocket  chamber.  One  technique  used  for  solid  propellants  is  the  magnetic 
flowmeter.  Micciand  Caveny^  determined  the  real  part  of  the  pressure-coupled  response  function 
in  a  forced  longitudinal  wave  motor  which  burned  solid  propellants.  A  spinning  gear  at  the  exit  of 
the  rocket  motor  nozzle  induced  acoustic  oscillations  (unsteady  pressure  and  velocity)  within  the 
motor.  The  experiments  proved  the  validity  of  using  a  magnetic  flowmeter  in  solid  rockets 
independent  of  the  propellant  type.  Wilson  and  Micci^  implemented  a  similar  experimental  setup  to 
measure  the  real  and  imaginary  parts  of  the  pressure-coupled  response  function.  Response 
functions  were  found  for  two  formulations  of  AP/HTPB  composite  propellants.  Cauty,  Comas, 
Vuillot  and  Micd'^  used  a  similar  experimental  setup  to  determine  the  real  part  of  the  pressure- 
coupled  response  function.  This  experiment  deduced  the  response  functions  from  the  steady  and 
unsteady  pressure  and  the  unsteady  velocity.  One  experiment,  completed  by  Holme^,  measured 
the  steady  velocity  in  the  exhaust  of  a  liquid  rocket  with  an  electromagnetic  flowmeter.  The 
flowmeter  used  the  alternating  magnetic  field  from  an  electromagnet.  Holme  concluded  that  the 
electromagnetic  flowmeter  worked  properly  for  several  types  of  liquid  propellants  excluding  a 
hydrogen  peroxide  monopropellant  due  to  its  low  flame  temperature. 

EXPERIMENTAL  APPROACH 


Sub-Scale  Rocket  Chamber 

A  sub-scale  rocket  chamber  previously  used  to  measure  liquid  oxygen  drop  size  and 
velocity®  was  modified  to  include  a  segment  which  could  measure  unsteady  pressure  and 
velocity^'®.  The  design  of  the  chamber  was  based  on  a  single  injection  element  of  the  Space  Shuttle 
Main  Engine  (SSME)  Prebumer.  Liquid  oxygen  was  injected  coaxially  inside  a  surrounding 
gaseous  hydrogen  stream.  The  chamber  was  modular  in  desi^  and  consisted  of  six  segments 
including,  the  injector,  ignitor  segment,  pressure  release  mechanism,  magnetic  flowmeter  segment, 
window  segment  and  nozzle.  The  total  length  of  the  chamber  was  0.377  m  (14.8  in)  and  the  width 
of  the  hexagonal  cross-section  from  one  face  to  the  other  was  0.028  m  ( 1 . 12  in). 

Acoustic  oscillations  were  induced  within  the  chamber  by  rotating  a  gear  inmiedi^ely 
downstream  of  the  nozzle  as  shown  in  Figure  1.  When  the  rotating  teeth  of  the  gear  partially 
obstructed  the  flow  from  the  converging  nozzle  the  throat  area  was  periodically  varied.  The  gear 
was  constructed  of  steel  and  contained  28  teeth.  The  height  and  width  of  the  teeth,  0.0127  m  (0.5 
in),  were  approximately  the  same  size  as  the  diameter  of  the  nozzle  throat,  0.0132  m  (0.52  in).  A 


2 


protective  steel  box  was  constructed  to  contain  and  divert  the  exhaust  gases  away  from  the  rocket 
chamber.  The  gear  was  driven  by  an  explosion-proof  motor  which  in  turn  was  driven  by  an 
invertor.  The  invertor  remotely  controlled  the  frequency  of  the  signal  transmitted  to  the  motor 
which  controlled  the  motor  speed. 

Magnetic  Flowmeter 

Unsteady  pressure  and  velocity  measurements  were  acquired  within  the  magnetic 
flowmeter  segment  of  the  rocket  chamber.  This  segment  contained  a  port  for  a  high  frequency 
pressure  transducer  to  measure  unsteady  pressure  and  a  magnetic  flowmeter  to  measure  unsteady 
velocity  at  0.139  m  (5.47  in)  downstream  of  the  injector.  Tbe  magnetic  flowmeter  consisted  of  a 
permanent  magnet  and  two  electrodes.  The  interaction  between  the  magnetic  field  and  the 
electrically  conducting  combustion  gas  within  the  chamber  induces  a  voltage  across  the  electrodes. 

From  electromagnetic  theory  the  electric  field,  £,  is  defined  as  E  =  UxB  where  U  is  the 
flow  velocity  of  an  electrically  conducting  fluid  and  B  is  the  magnetic  field.  When  two  metal  rods 
are  placed  in  the  electric  field  perpendicular  to  the  direction  of  a  gas  flow  a  potential  difference  is 
generated.  This  potential  difference  is  proportional  to  the  product  of  the  flow  velocity  and  the 
magnetic  field  according  to  Faraday’s  Law  .  V  =  OJLB .  In  this  equation  is  a  proportionality 
constant  of  the  order  1 ,  L  is  the  distance  between  the  electrodes  and  V  is  the  induced  voltage. 

A  cross  sectional  view  of  the  magnetic  flowmeter  used  in  the  experiments  is  shown  in 
Figure  2.  The  path  of  the  combusting  gases  passed  through  the  magnetic  flowmeter  segment 
perpendicular  to  the  page.  The  magnetic  field  was  produced  by  a  horseshoe  permanent  magnet 
with  a  gap  density  of  0.19  Tesla.  The  electrodes  protruded  into  the  chamber  approximately  0.003 
m  (0. 12  in)  with  the  distance  between  them  approximately  0.022  m  (0.88  in). 

Experimental  Procedure 

The  hot  fire  combustion  tests  were  conducted  at  the  Cryogenic  Combustion  Laboratory  at 
the  Pennsylvania  State  University.  For  each  test,  high  frequency  pressure  and  velocity  data  were 
acquired  from  the  piezoelectric  pressure  transducer  and  magnetic  flowmeter.  The  gear  which 
modulated  the  exhaust  was  rotated  at  a  frequency  between  2500  and  2700  teeth  per  second.  These 
frequencies  were  near  the  first  longitudinal  resonance  frequency  of  the  chamber.  The  mean 
pressure  for  the  tests  ranged  from  2.6  to  3.0  MPa  (375  to  425  psi).  Table  1  lists  the  actual 
operating  conditions  of  each  test. 

Table  1  -  Operating  conditions  of  the  experiments 


Test  No. 

Pressure 

Mixture  VelocityMomentum 

Liquid 

Weber  No.  Frequency 

(psi) 

(MPa) 

Ratio 

Ratio 

Ratio 

Reynolds  No 

• 

(Hz) 

Test  1 

379 

2.61 

4.08 

19.34 

0.91 

788,556 

441,228 

2,669 

Test  2 

400 

2.76 

3.82 

20.16 

1.01 

774,704 

463,969 

2,540 

Tests 

414 

2.86 

3,87 

19.15 

0,95 

738324 

377,074 

2,613 

Test  4 

429 

2.96 

4.94 

14.24 

0.55 

830.486 

271394 

2,607 

Test  5 

405 

2.80 

5.50 

14.04 

0.49 

856305 

252300 

2,624 

Test  6 

410 

2.83 

5.29 

14.50 

0.53 

899.763 

316,532 

2,568 

EXPERIMENTAL  RESULTS 


The  voltage  difference  measured  by  the  magnetic  flowmeter  was  converted  into  velocity  as 
described  in  the  previous  section.  The  dimensioidess  proportionality  constant,  a,  could  not  be 
determined  and  was  assumed  to  be  0.5.  For  most  of  the  tests  the  pressure  signal  oscillated 
between  ±0.28  MPa  (±40  psi)  and  the  velocity  signal  had  an  average  value  between  ±1.5  m/s 
(dt5  ft/s).  The  pressure  and  velocity  time  records  were  transformed  into  the  frequency  domain  by  a 
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Fast  Fourier  Transform  (FFTy®  to  calculate  the  Power  Spectrum  Density  (PSD)  of  the  pressure, 
,  and  velocity,  S„(f) ,  signals  in  the  frequency  domain.  Hgures  3  and  4  are  plots  of  the  PSD 
of  the  pressure  and  velocity  signals  for  a  test  with  a  constant  gear  frequency  of  2669  teeth  per 
second.  The  data  is  plotted  in  three-dimensions  on  a  waterfall  plot  to  show  dependence  with  time. 
Each  horizontal  line  on  the  plot  represents  data  taken  over  0. 178  seconds.  The  pressure  PSD 
clearly  shows  the  acoustic  oscillations  in  the  chamber  at  the  frequency  induced  by  the  gear.  Both 
the  first  and  second  harmonics  are  visible  in  the  graph  at  approximately  2670  and  5340  Hz.  As 
expected  the  velocity  PSD  is  much  noisier,  however  the  induced  frequency  is  stiU  visible  at  both 
the  first  and  second  harmonics  corresponding  to  the  pressure  wave. 

The  pressure  and  velocity  signals  were  correlated  to  obtain  the  Cross  Spectrum  Density 
(CSD)  between  the  pressure  and  velocity  signals  in  the  frequency  domain.  Figure  5  is  a  wateifaU 
plot  of  the  CSD  for  the  same  test  shown  in  Figures  3  and  4.  The  CSD  shows  the  pressure  and 
velocity  signals  to  be  related  at  the  gear  forcing  frequency  and  the  plot  is  much  less  noisy  than  the 
waterfall  plot  of  the  velocity  PSD.  The  CSD  provided  a  means  to  calculate  the  phase  difference 
between  the  pressure  and  velocity  signals.  The  phase  was  calculated  from  the  real  and  imaginary 
parts  of  the  complex  CSD  frequency  signal 

The  error  associated  with  performing  the  CSD  is  determined  by  the  coherence  function  as 
defined  in  Equation  1  where  ft  is  the  frequency  corresponding  to  the  gear  forcing  frequency.  The 
coherence  is  a  ratio  of  the  averaged  CSD  to  the  product  of  the  measured  pressure  PSD  and  the 
measured  velocity  PSD.  The  coherence  lies  between  the  values  of  zero  and  one  and  equals  one  for 
a  noiseless  signal.  Measured  and  averaged  data  have  a  coherence  less  than  one.  The  coherence 
predicts  the  error  in  the  CSD  phase  calculation.  Equation  2  defines  the  phase  error  where  n,  is  the 
number  of  points  over  which  the  phase  is  averaged. 


ACOIJSTICMODEL 


(1) 

(2) 


Governing  Equations 

A  linear  one-dimensional  acoustic  model  was  derived  to  numerically  solve  for  the  mean  and 
unsteady  pressure  and  velocity  within  a  distributed  combustion  rocket  chamber  following  the  work 
of  Med'  ^  for  a  solid  rocket  motor.  The  model  starts  with  the  governing  equations  for  a  one¬ 
dimensional,  unsteady,  compressible  flow  in  a  constant  area  rocket  chamber.  These  equations 
include  the  continuity  equation,  conservation  of  momentum  and  energy,  and  the  perfect  gas  law 
with  mass  and  energy  addition  to  describe  the  fluid  properties  of  the  products  of  combustion.  The 
mass  and  energy  addition  is  described  by  a  mass  flux  term  and  is  determined  by  the  amount  of 
liquid  oxygen  injected  into  the  rocket  chamber. 

The  governing  equations  are  linearized  by  assuming  each  fluid  property  is  the  sum  of  a 
mean  and  a  small  perturbation  value  and  is  assumed  to  be  sinusoidal  in  time.  Terms  of  the  order 
M  or  less  are  neglected.  The  unsteady  adiabatic  flame  temperature,  fy ,  is  described  as  the  sum  of 
an  isentropic  and  non-isentropic  unsteady  temperature  contribution  as  shown  in  Equation  3 .  The 
pressure-  and  velocity-coupled  response  functions  for  the  liquid  propellant  evaporation  and 
combustion  are  given  in  Equation  4. 
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(3) 

(4) 


mi,  Tf 


Solution  Technique 

The  five  linearized  equations  given  by  the  derivation  of  the  model  were  solved 
simultaneously  for  different  pressure-  and  velocity-coupled  response  functions  given  initial 
conditions.  These  initial  conations  were  set  for  the  properties  at  the  chamber  wall  near  the 
injector.  The  initial  condition  for  the  mean  density  equaled  the  density  of  the  hydrogen  entering  the 
chamber  from  the  injector.  The  mean  pressure  was  set  to  the  mean  chamber  pressure.  The 
unsteady  velocity  equaled  zero  at  the  injection  wall  of  the  chamber.  The  mean  velocity  was  based 
on  the  injection  velocity  of  the  hydrogen  into  the  chamber.  The  unsteady  pressure  was  assumed  to 
be  approximately  10%  of  the  mean  chamber  pressure.  These  initial  conditions  are  summarized  in 
Table  2. 


In  order  to  predict  the  conditions  in  the  chamber  the  reaction  was  simulated  with 
CHEMKIN'^  a  well-stirred  reactor  model.  The  inputs  included  the  chamber  pressure,  total 
propellant  mass  flow  rate,  chamber  volume,  fuel  mass  fraction,  oxidizer  mass  fraction  and 
equivalence  ratio.  The  program  predicted  the  adiabatic  flame  temperature  and  the  mass  fraction  of 
the  product  species  which  allowed  for  the  calculation  of  the  gas  properties,  including  the  mixture 
molecularweight,  the  mixture  gas  constant,  the  ratio  of  specific  heats  and  the  speed  of  sound  in  the 
chamber. 


The  remaining  parameters  needed  to  solve  the  differential  equations  were  the  mass  flux, 
frequency,  and  pressure-  and  velocity-coupled  response  functions.  The  mass  flux  source  term  in 
the  governing  equations  was  defined  as  the  mass  flow  rate  of  LOX  per  chamber  volume.  The 
mass  flux  was  assumed  to  be  a  parabolic  function  of  downstream  distance.  The  frequency  was 
defined  from  the  frequency  corresponding  to  the  maximum  power  on  the  CSD  waterfall  plots 
which  was  also  the  forcing  frequency  of  the  gear.  The  pressure-  and  velocity-coupled  response 
functions  were  varied  and  will  be  discussed  more  in  the  next  section. 


RESULTS  FROM  MODEL 


The  model  was  solved  for  various  real  values  of  the  pressure-coupled  response  function, 
Rp ,  and  the  velocity  coupled  response  function,  R, .  The  results  from  the  model  predicted  the 
phase  difference  between  the  unsteady  pressure  and  velocity.  The  model  was  related  to  the 
experiment  by  comparing  the  phase  difference  of  the  model  at  0.139  in  (5.47  in)  downstrwm  of 
the  injector  face  to  the  phase  difference  calculated  from  the  CSD  for  the  experiments.  This  distance 
of  0.14  m  corresponds  to  the  downstream  position  of  the  pressure  transducer  and  magnetic 
flowmeter  in  the  rocket  chamber. 

Table  2  -  Initial  conditions  for  linear  acoustic  model 


Propert 

y 

Initial 

Guess 

p 

p 

p^ 

u 

0 

u 

p 

O.IOF^ 
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A  parametric  study  was  completed  by  varying  and  R,  for  each  test.  Real  integer  values 
were  used  for  each  response  function:  R^  was  varied  between  0  and  5  and  R,  ranged  from  0  to  1 2. 
As  the  value  of  the  response  functions  changes,  the  downstream  location  where  the  phase 
difference  shifts  from  90®  to  -90°  also  changes  producing  a  range  of  values  for  the  phase  difference 
at  0. 14  m.  The  phase  difference  shifts  from  a  positive  value  at  R  =  0  to  a  negative  value  at  R  =  12 . 
For  values  of  R^,  near  zero  the  phase  difference  shifts  from  positive  to  negative  through  ±180° 
whereas  for  values  of  near  five  the  phase  shifts  through  0°.  This  particular  phenomenon  is 
observed  because  the  model  predicts  the  node  of  the  pressure  wave  and  hence  the  maximum 
unsteady  velocity  magnitude  to  be  near  the  position  of  the  magnetic  flowmeter.  By  changing  the 
values  of  the  pressure-  and  velocity-coupled  response  functions  the  pressure  node  and  velocity 
maximum  may  be  moved  slightly.  This  slight  movement  can  strongly  affect  the  phase  difference 
between  the  unsteady  pressure  and  velocity. 

The  experimental  phase  difference  for  each  test  was  chosen  based  upon  its  coherence. 
Points  corresponding  to  the  highest  coherence  or  the  least  amount  of  statistical  error  were 
compared  to  the  phase  difference  computed  by  the  model.  By  comparing  the  phase  differences  for 
each  test  to  the  model  a  range  of  values  for  R^  and  R,  were  found.  For  a  phase  difference  of  -60°, 
the  model  predicts  a  match  where  =  3,  4  or  5  and  R,  is  between  4  and  8 . 

The  range  for  the  pressure-coupled  response  function  for  each  test  is  plotted  in  Figure  6. 
Test  1  in  Figure  6  corresponds  to  the  test  discussed  in  the  above  text.  The  conditions  of  Tests  1,  2 
and  3  are  approximately  the  same  as  shown  previously  in  Table  1 .  Likewise  Test  4,  5  and  6  were 
run  at  similar  conditions.  The  pressure-coupled  response  function  for  the  conditions  of  Tests  4,  5 
and  6  is  between  the  range  of  zero  and  two.  Although  the  range  for  R^,  is  quite  large  at  any  given 
frequency,  for  Tests  1,  2  and  3,  R^  changes  from  between  zero  and  two  at  2540  Hz  to  between 
three  and  five  at  2670  Hz. 

A  similar  plot  for  the  velocity-coupled  response  function  is  shown  in  Hgure  7.  This  graph 
indicates  that  R  decreases  from  between  two  and  four  at  2540  Hz  to  zero  from  2610  to  2670  Hz 
for  Tests  4,  5  and  6.  The  velocity-coupled  response  function  ranges  from  four  to  eight  for  the 
frequency  range  between  2540  and  2670  Hz  for  T ests  1 ,  2  and  3 . 

CONCLUSIONS 

A  rotating  gear  near  the  exit  of  the  nozzle  of  a  sub-scale  liquid  rocket  chamber  induced 
acoustic  oscillations  within  the  chamber.  Pressure  transducers  measured  the  mean  and  unsteady 
pressures  within  the  chamber.  A  magnetic  flowmeter  equipped  with  a  permanent  magnet  measured 
the  unsteady  velocity  of  the  combusting  gases.  The  experiments  have  proven  that  the  magnetic 
flowmeter  is  capable  of  acquiring  unsteady  velocity  data  within  a  gaseous  hydrogen/liquid  oxygen 
rocket  chamber. 

The  unsteady  pressure  and  velocity  time  signals  were  transformed  into  the  frequency 
domain  with  a  Fast  Fourier  Transform.  The  two  frequency  signals  were  then  correlated  to 
determine  the  frequency  at  which  both  of  the  signals  were  oscillating,  i.e.  the  forcing  frequency  of 
the  rotating  gear  in  the  experiment.  After  averaging  the  phase  difference  between  the  unsteady 
pressure  and  velocity  was  determined  at  this  frequency.  TTie  phase  error  was  calculated  from  the 
coherence  function  at  the  forcing  frequency. 

A  linearized  one-dimensional  acoustic  model  was  developed  to  deduce  the  pressure-  and 
velocity-coupled  response  functions  corresponding  to  the  experimental  results.  The  equations 
derived  from  the  model  were  solved  for  the  mean  density,  pressure  and  velocity  and  the  unsteady 
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pressure  and  velocity..  The  flow  properties  measured  during  the  experiment  were  input  into  the 
model  as  initial  conditions. 

The  model  predicted  the  phase  difference  between  unsteady  pressure  and  velocity  for 
various  real  valued  pressure-  and  velocity-coupled  response  functions.  Response  functions  were 
found  where  the  phase  difference  predicted  by  the  model  matched  the  phase  difference  calculated 
from  the  experiment.  Each  experimental  phase  difference  chosen  produced  a  range  of  values  for 
the  pressure-  and  velocity-  coupled  response  functions.  The  real  part  of  pressure-coupled  response 
fimction  ranged  from  zero  to  five  and  the  real  part  of  the  velocity-coupled  response  fimction  ranged 
from  zero  to  twelve. 
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Figure  1  -  Gear  and  nozzle  setup 


Figure  2  -  Magnetic  flowmeter  chamber  segment 


Figure  3  -  Unsteady  pressure  power  spectral  density 


Figure  4  -  Unsteady  velocity  power  spectral  density 


Hgure  5  -  Unsteady  pressure  and 
velocity  cross  spectral  density 


Rgure  6  Pressure-coupled  response  function  at  different 
frequencies  and  test  conditions. 


Rgure  7  Velocity-coupled  response  function  at  different 

frequencies  and  test  conditions. 
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